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Abstract 

It has been shown that in the twisted mass formulation of Lattice QCD 
at maximal twist large cutoff effects are generated when the quark mass 
becomes of 0{aK^ CD ). In general, these effects can be suppressed in two 
ways: either by choosing the critical quark mass in an "optimal way" , or by 
adding the Clover term to the twisted action. We investigate the second op- 
tion by performing a quenched lattice QCD simulation with twisted Clover 
fermions and pion masses as low as 280 MeV. We show that the Clover term 
is indeed efficient in reducing the large cutoff effects. In particular, the so- 
called bending phenomenon observed in the determination of the pion decay 
constant is cured in this way. In addition, by using maximally twisted Clover 
fermions, we provide a non-perturbative determination of the vector current 
renormalization constant Zy as well as of the non-perturbatively renormal- 
ized light quark masses. Finally, we calculate the connected contribution 
to the charged-neutral pseudoscalar meson mass splitting, finding that the 
introduction of the Clover term in the twisted action is also beneficial, in 
the quenched approximation, in reducing cutoff effects related to the isospin 
symmetry breaking at finite lattice spacing. 



1 Introduction 



The twisted mass formulation of Lattice QCD (tmLQCD) [1, 2] offers a number of 
interesting advantages. The absence of real unphysical zero modes, which affect the 
standard Wilson discretization of the fermionic action, prevents the appearance 
in this framework of the exceptional configurations in quenched and partially- 
quenched numerical simulations. This allows to investigate QCD dynamics with 
much lighter valence quarks. In addition, twisting the Wilson term in the action (in 
the physical basis) considerably simplifies the renormalization pattern of several 
local operators. Notable examples in this respect are the determination of the pion 
decay constant, which can be performed without introducing any renormalization 
factor [1], and of the kaon Bk parameter, which does not require subtraction 
of wrong chirality operators [L3]. 1 Furthermore, tmLQCD with a maximal twist 
angle (a = vr/2), denoted as maximally twisted mass LQCD (Mtm-LQCD), ensures 
automatic (or almost automatic) (9(a)-improvement for the physically interesting 
quantities computed on the lattice [6]. 

The drawback, however, is that with tmLQCD both parity and flavor symme- 
tries are explicitly broken at finite lattice spacing. The breaking of parity generates 
a mixing with wrong parity states which may affect, in some relevant cases, the 
large time behavior of lattice correlation functions. The breaking of flavor sym- 
metry is responsible for generating a mass splitting between otherwise degenerate 
isospin hadronic partners. Although the breaking effects are parametrically of 
O (a 2 ) at maximal twist, this splitting has been found to be numerically significant 
in the case of the neutral and charged pseudoscalar mesons [7]. 

Another important issue of practical interest for Mtm-LQCD is the presence 
of large cutoff effects which are generated when the quark mass m q becomes of 
0(ciAq CD ) or smaller. This phenomenon has been the subject of recent numerical 
and theoretical investigations. From the point of view of the Symanzik expansion, 
these large cutoff effects have been studied in Ref. [8]. In correlation functions 
they are related to multiple insertions of the dimension-five operator £ 5 entering 
the Symanzik effective Lagrangian for Mtm-LQCD close to the continuum limit. 
Since, at maximal twist, £5 has the quantum numbers of the neutral pion, its 
insertions can create from the vacuum Goldstone bosons with zero momentum, 
thus generating discretization effects which are proportional, at the leading order, 
to (a/M 2 ) 2 ~ (a/m q ) 2 . In this situation, the value of the quark mass is limited to 
the region m q ^ ah 2 Q CD . Being proportional to (a/m q ) 2 , such lattice artifacts have 
been also called "infra-red divergent" (IRD) discretization effects. Note, however, 
that they actually do not represent any infra-red divergence, since the chiral limit 

1 For an alternative approach which avoids this subtraction even when working with ordinary 
Wilson fcrmions see Refs. [4,5]. 
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can be taken only after the continuum extrapolation has been performed. 

A further step forward in understanding the IRD cutoff effects has been made in 
Ref . [9] , where it is shown that their contribution to the Symanzik expansion can be 
explicitly resummed using twisted-mass chiral perturbation theory (tmxPT) [10]. 
For a simple quantity like the pion decay constant f n , the leading order result of 
tm^PT reads 



where / is the pion decay constant in the chiral limit and fx and 8m represent 
the (renormalized) twisted and non-twisted quark masses respectively. Since at 
maximal twist 8m is a quantity of 0(a) introduced by a discretization error in the 
determination of the critical mass, Eq. (1) shows that the deviation of f n from its 
value in the chiral limit has the form of an IRD cutoff effect. The leading correction 
is proportional to (8m/ fx) 2 ~ (a/M 2 ) 2 and its contribution becomes sizable when 
the twisted mass fx is small and comparable to Sm, i.e., fx ~ 0(aAn CD ). 

The requirement fx ^ ak 2 Q CD is in practice rather restrictive. For a typical 
inverse lattice spacing a" 1 ~ 2 GeV, it implies pion masses M w > 500 MeV. In 
order to allow for simulations with smaller values of the quark masses (down to 
fx <; a 2 Is?Q CD ), and keep discretization effects under control, two strategies have 
been proposed [8] . Either one uses an "optimal" determination of the critical mass 
m cr or the (twisted) Clover term is introduced in the quark action. In the latter 
case, the coefficient csw has to be fixed to the same value that guarantees on-shell 
(9(a)-improvement with ordinary Wilson fermions. 

The idea behind the first strategy, namely the optimal determination of the 
critical mass is to tune the 0(a) contribution to m cr in such a way that the 
dimension-five operator in the Symanzik effective Lagrangian has vanishing matrix 
element between the vacuum and the single pion state, i.e. (0 1 i2s 1 7r°) = [8]. 
Since the IRD cutoff effects have been shown to be proportional at the leading 
order to the previous matrix element, this procedure guarantees their suppression 
and allows much smaller values of the quark masses to be safely reached in the 
simulation. In practice, this procedure corresponds to fix the value of the critical 
mass by requiring the vanishing of the PCAC quark mass in the twisted basis, as 
it was anticipated in Refs. [11, 12]. In the following, we will denote as m°^ this 
"optimal" determination of the critical mass. 

The numerical studies reported in Refs. [13, 14] show that such a strategy is 
indeed efficient. In particular, in agreement with theoretical expectations, it helps 
in suppressing the pronounced bending behaviour of the pion decay constant, i.e., 
the appearance of large lattice artifacts in this quantity when the quark mass be- 
comes of 0(<xA?q CD ) or smaller. Note, however, that one has to make a preliminary 




(1) 
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numerical study, with twisted fermions, to determine the optimal value of the crit- 
ical mass. Such a determination may be computationally expensive, particularly 
in the case of unquenched simulations. 

In this paper we explore the second strategy of Ref. [8], namely the simulation 
of Mtm-LQCD with the inclusion of the Clover term in the action. In this case, 
in the Symanzik expansion of the lattice effective Lagrangian at 0(a) only a term 
proportional to the square of the twisted quark mass survives [15], that is £5 = 
0{ajj?) = 0(aM%). Thus, the introduction of the Clover term is expected to be 
efficient in cancelling the leading order IRD cutoff effects as much as the tuning of 
the critical mass in the twisted Wilson case. 

From the previous discussion it should be clear that the Mtm-Clover approach 
does not require any specific prescription for the determination of the critical mass. 
We show here that this is indeed the case. We have performed two numerical sim- 
ulations with Mtm-Clover fermions and with a critical quark mass determined in 
the non-twisted theory either by requiring the vanishing of the pion mass in the 
chiral limit (denoted in the following as m^ on ) or by imposing the vanishing of the 
PCAC mass (m^ r CAC ). In both cases, we do not observe any bending phenomenon 
in the pion decay constant, even at values of quark masses smaller than 0(aAq CD ). 
We also compare our results for the pion mass and decay constant with those ob- 
tained by the %LF Collaboration [13] using Mtm- Wilson fermions with the optimal 
choice of the critical mass, finding a very good agreement. The appealing feature 
of the Mtm-Clover approach is that the determination of the critical mass does not 
require any specific tuning and can be obtained from a single simulation performed 
with non-twisted Clover fermions. 

Besides investigating the effectiveness of Mtm-Clover fermions in suppressing 
the IRD cutoff effects, we also present a non-perturbative determination of the local 
vector current renormalization constant Zy which renormalizes, in Mtm-LQCD, 
the local axial- vector current (in the physical basis). This constant is expected 
to be identical, up to 0(a 2 ) effects, to the value determined by using non- twisted 
Clover fermions. Indeed, we find a very good agreement between our result and the 
non-twisted determination of Zy performed in Ref. [16] with the non-perturbative 
RI-MOM method. 

In addition, we present a calculation of the strange and of the average up/down 
quark masses using the value of the quark mass renormalization constant deter- 
mined non-perturbatively in Ref. [16]. With respect to previous quenched determi- 
nations of light quark masses, here we benefit from the use of Mtm-LQCD with the 
Clover term to probe ever lower values of the quark masses, and thereby reduce 
the errors due to the chiral extrapolation when determining the light up/down 
quark masses on the lattice. 

Finally, we investigate the isospin symmetry breaking effects, induced at finite 
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lattice spacing by the tmLQCD action, by studying the contribution of connected 
diagrams to the charged-neutral pion mass splitting. This connected contribution 
can be interpreted as the charged-neutral mass splitting between pseudoscalar 
mesons composed by quarks and antiquarks belonging to different, though degen- 
erate, isospin doublets. By comparing the results obtained with Mtm- Wilson and 
Mtm-Clover fermions, we find that the introduction of the Clover term is beneficial 
also in reducing these cutoff effects, at least for the pure connected contribution 
considered in this study and in the quenched approximation. 

It is also worth noting that the same 0(a 2 ) operator that provides the leading 
contribution to the pion mass splitting in Mtm-LQCD also affects the determina- 
tion of both the charged and the neutral pion masses when working with ordinary 
Wilson/Clover fermions. In this respect, therefore, the use of Mtm-LQCD does 
not introduce any additional uncertainty. 

The plan of the paper is as follows. In Section 2 we give the details of the lattice 
simulation and present our results for the pseudoscalar meson mass and decay 
constant; we show the effectiveness of the introduction of the Clover term in curing 
the bending phenomenon. We present the determination of the renormalization 
constant Zy in Section 3.1 and the calculation of the light quark masses in Section 
3.2. In Section 3.3 we evaluate the contribution of connected diagrams to the 
charged-neutral pion mass splitting and compare the results with those obtained 
by using Mtm- Wilson fermions with the optimal choice of the critical mass. Our 
conclusions are summarized in Section 4. 



2 Pseudoscalar meson masses and decay constants 
with the Mtm-Clover action 

The numerical results presented in this study have been obtained from a set of 300 
gauge configurations generated, in the quenched approximation, with the standard 
Wilson gauge action at f3 — 6.0, which corresponds to an inverse lattice spacing 
a 1 ~ 2.0 GeV. The lattice volume is I6 3 x 32 so that the physical length of the 
lattice is about 1.6 fm in the spatial directions. Using this set of gauge configu- 
rations, quark propagators have been computed by implementing three different 
choices for the quark action: 

a) A simulation with Mtm-Clover fermions and the critical mass determined 
from the PCAC relation in the non-twisted theory, m cr = m^ r CAC . The 
value of the critical hopping parameter, fi^ r CAC = 0.135217, corresponding 
to flm^. CAC ~ —0.3022, has been obtained by using the results of a previous 
simulation done by our collaboration with standard Clover fermions at the 
same value of the gauge coupling on the volume 24 3 x 64 [16]. 
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b) A simulation with Mtm-Clover fermions and the critical mass determined 
from the vanishing of the pion mass in the non-twisted theory, m cr = m^ on . 
The value of the critical hopping parameter, n^° n = 0.135293, corresponds to 
am^ on ~ —0.3043 and has been determined from the same set of non-twisted 
Clover data used to determine the value of m^, CAC in simulation a). 

c) A simulation with Mtm- Wilson fermions with the critical mass determined 
in the optimal way, that is m cr = m° Pt . We used the value = 0.157409 
obtained by the xLF Collaboration from a simulation at f3 — 6.0 on the 
volume 16 3 x 32 [13]. 

In simulations a) and b) with Mtm-Clover fermions the coefficient of the Clover 
term has been fixed to the value determined non-perturbatively in Ref. [17], namely 
csw = 1-769. 

In all three sets of simulations we have computed quark propagators with the 
same 9 values of the bare twisted quark mass a/i (cf. Table 1), which are equal to 
those used in Ref. [13]. The range of chosen a/i values covers pion masses between 
approximately 280 MeV and 1.1 GeV. In order to invert the Dirac operator for such 
a broad set of quark masses we employed the multiple mass solver algorithm [18]. 

To examine the impact of cutoff effects when working with maximally twisted 
fermions at small values of the quark mass, we investigate whether or not the 
inclusion of the Clover term is efficient in curing the bending phenomenon observed 
in the determination of the pion decay constant with Mtm- Wilson fermions, when 
the critical mass is determined in a non-optimal way. 

Pseudoscalar meson masses and decay constants have been evaluated by study- 
ing the large time behavior of the following two-point correlation functions: 

C PP (t) = J2(P^^)Pk0)) , C AP (t) = J2(Mt,x)PK0)) , (2) 

X X 

where P 5 = u'-f^d and = u^^^d. Throughout this paper we choose to work in 
the physical basis and use mass-degenerate u and d quarks. 

Using the completeness relation and considering large enough time separations, 
one gets 

Cpp[t) «>l if e ~ MpT ' 2 cosh (Mp(r/2 " ()) • 

f>t ^^e-^ S inh(Mp(T/2-t)) , (3) 

where M P is the mass of the lightest pseudoscalar meson state, \JZ PS = | (0|P 5 (0) |vr) | 
and f P M P = Z v (0\A (0)\ir). Z v is the renormalization constant of the local 
axial-vector current with Mtm-LQCD, which corresponds to the renormalization 
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constant of the local vector current with the non-twisted quark action [6]. The 
pseudoscalar decay constant can therefore be extracted by studying the ratio 
CAp(t)/Cpp(t) at large time distances, 



, rv V Zps 



^§-coth(M P (T/2-t)) 



(4) 



*>i 



h = 2^ ■ (6) 



An alternative determination of the decay constant [1] , which does not require 
the introduction of any renormalization constant, can be obtained by exploiting 
the consequences of the axial-vector Ward identity (AWI) for Mtm-LQCD. Up to 
O (a 2 ) correction, the relevant identity reads 

Z v J2(doA (t,x)4(0)) = 2/x ^(PsM^O)) , (5) 

X X 

where \i is the bare twisted mass. Combining the above relation with Eqs. (3) one 
gets 

V Zps 

where both Zps and Mp can be determined from the correlation function Cpp(t) 
of Eq. (3). We have evaluated the pseudoscalar meson decay constant fp using 
both the "direct" method (4) and the "indirect" one (6), obtaining always a very 
good agreement within the statistical errors. 2 For the direct method, we used 
the non-perturbative RI-MOM determination z^ lover ^ = 0.772 [16] and the one- 
loop boosted perturbative estimate z^ dson ^ = 0.63 [19,20] in the Mtm-Clover 
and Wilson case respectively. In what follows, we only quote the results for fp 
obtained from the indirect method (6). 

Our result for the pseudoscalar meson masses and decay constants obtained 
from simulations a), b) and c) are collected in Table 1. These results are also 
shown in Figs. 1 and 2, as a function of the bare quark mass a fx. The Mtm- Wilson 
results from our simulation c) agree well with those presented in Ref. [13]. For 
comparison, we also show in these plots the results obtained in [13] by employing 
Mtm- Wilson fermions with the non-optimal choice of the critical mass, as well 
as the results obtained with non-twisted Clover fermions (at j3 = 6.0) using data 
produced by our collaboration for the studies described in Refs. [16,21]. In the 
non-twisted case, the quark mass a/i in Figs. 1 and 2 represents the bare AWI 
quark mass am AWI [21] multiplied by the axial- vector renormalization constant 

2 This agreement is further improved if the factor Mp in the denominator of Eq. (6) is replaced 
with Mp sinh(aMp)/a, where the sinh is originated from the discretized symmetric lattice version 
of the derivative in Eq. (5). Our results from Eq. (6) have been obtained with this replacement 
always implemented. 
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a \x 



Fig. 2: Same as in Fig. 1, but for the pseudo-scalar decay constant. The lines represent 
quadratic fit in a[i except for the results obtained with Mtm-Wilson fermions and the 
mcr° n determination of the critical mass. In this case, the fit also takes into account the 
tmxPT prediction expressed by Eq. (1 ). 
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a) Mtm-Clover, m^, CAC 


b) Mtm-Clover, m£ on 


c) Mtm- Wilson, m°f 


a/i 


alVlp 


ajp 


alVlp 


ajp 


alVlp 


ajp 


0.0038 


0.1359 (66) 


0.0690 (30) 


0.1308 (70) 


0.0739 (34) 


0.1233 (68) 


0.0684 (29) 


0.0057 


0.1607 (52) 


0.0710 (22) 


0.1571 (51) 


0.0741 (23) 


0.1515 (52) 


0.0696 (18) 


0.0076 


n 1 81 1 ( AR\ 
U.loll \fiX>) 


n 0797 

U.U I LI (lo) 


n 1 787 (A 7\ 
U.l lot [41 ) 


n n7/i 8 (~\q\ 
u.U(4o 


n 1 7^1 (A A \ 

U.l (01 [fi^) 


n n7i n (~\ a \ 

U.U f 1U ^ 141 J 


0.0113 


0.2150 (39) 


0.0754 (15) 


0.2137 (39) 


0.0767 (15) 


0.2127 (36) 


0.0739 (11) 


U.U101 


0.2449 (35) 


0.0779 (13) 


0.2444 (35) 


0.0787 (13) 


0.2445 (31) 


0.0767 (10) 


0.0302 


0.3401 (26) 


0.0865 (11) 


0.3404 (26) 


0.0868 (11) 


0.3425 (22) 


0.0861 (09) 


0.0454 


0.4182 (22) 


0.0942 (11) 


0.4187 (22) 


0.0944 (11) 


0.4223 (18) 


0.0944 (08) 


0.0605 


0.4873 (20) 


0.1012 (10) 


0.4879 (20) 


0.1014 (10) 


0.4930 (16) 


0.1019 (08) 


0.0756 


0.5512 (17) 


0.1077 (10) 


0.5519 (18) 


0.1079 (10) 


0.5585 (14) 


0.1092 (08) 



Table 1: Values of the bare twisted quark mass afi, of the pseudoscalar meson mass aMp 
and of the pseudoscalar meson decay constant afp from the three sets of simulations 
performed in this study. The quoted errors are purely statistical. 



Za- Indeed, this is the quantity that renormalizes with the same renormalization 
constant = l/Z P of the bare twisted mass a/i in the twisted formulation of 
LQCD. 3 

From the plots in Figs. 1 and 2 we observe that: 

• There is a very good agreement among the results for the pseudoscalar me- 
son masses calculated using various quark actions, all of which are 0(a)- 
improved. This may signal the smallness of discretization effects on this 
quantity at the value of the lattice spacing considered in the present study. 
Further investigations at different lattice spacings could better clarify this 
point. 

• IRD cutoff effects are clearly visible in the pseudoscalar meson decay constant 
at low quark masses (a/i ^ 0.02 ~ a 2 Ag CD ) when computed by using Mtm- 
Wilson fermions with the m^ on determination of the critical mass. This is 
the bending phenomenon observed in Ref. [13]. As shown in Fig. 2, these 
effects are strongly reduced either by choosing the optimal value of m cr , or 
by introducing the Clover term in the quark action. The reduction of IRD 
effects is almost the same within the two approaches. This is the main result 
of the present study. 

3 It should be also mentioned that the bare twisted quark mass an with Mtm-Clover fermions 
renormalizes differently from its counterpart with Mtm- Wilson fermions. However, the ratio of 
the corresponding renormalization constants, ^( 3 c(ol ' er )^( 3 H/dso ") j j s ver y c l ose to unity (it is 
0.97 when using the one- loop boosted perturbative expressions). 
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• As shown in Fig. 1, IRD cutoff effects are negligible in the case of the pseu- 
doscalar meson mass. This is consistent with the observation that IRD cutoff 
effects in the determination of this mass are softened by an additional factor 
of M 2 P [8]. 

• The results for the pseudoscalar meson masses obtained with the Mtm-Clover 
action and the m*£ on determination of the critical mass are not shown in 
Fig. 1, because they are practically indistinguishable from those obtained 
with m^ r CAC (cf. Table 1). In the case of the decay constant, instead, the 
two Mtm-Clover determinations with different choices of m cr show a small 
spread at very low quark masses. This may signal the contribution of higher- 
order IRD cutoff effects that we cannot cure. 

For illustrative purposes we quote the results for the pion decay constants 
obtained by quadratically extrapolating the lattice data to the physical pion mass. 
We obtain f w = {138(8), 144(9), 138(8)} MeV in simulations a), b) and c) 
respectively, to be compared with the experimental value f w = 132 MeV. 



3 More numerical results 



3.1 The renormalization constant Zy 

A non-perturbative determination of the renormalization constant Zy 4 can be 
obtained using the two-point AWI of Eq. (5) as 

7 -i _ (C AP (t + l)-C AP (t-l))/2 

Zv ~ 2f,C PP (t) ' [n 

where in the numerator we used the symmetric (O(a)-improved) version of the 
lattice time derivative. 

We have computed Zy from Eq. (7) using both the Wilson and the Clover 
action at maximal twist. In the Clover case, the results obtained from simulations 
a) and b), which only differ by the choice of the critical mass, are in good agreement 
within the statistical errors. For this reason, only data from simulation a) and c) 
will be discussed below. Our results for Z v are shown in Fig. 3 as a function of 
the bare twisted quark mass. 

The dependence of Zy on the quark mass is induced by discretization effects, 
which are predicted to be of 0(a 2 ) at maximal twist. Thus one can expect both 

4 In this paper we follow the standard practice to denote with Zy the renormalization constant 
of the axial- vector current in Mtm-LQCD. The reason is that the axial- vector current corresponds 
to the vector current in the twisted basis, which also implies that Zy coincides with the vector 
current renormalization constant of the standard non-twisted theory. 
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0.60 



0.000 



0.002 0.004 
(a n) 2 



0.006 




0.08 



Fig. 3: Values of the renormalization constant Zy as a function of the bare twisted quark 
mass as obtained by using the Mtm-Wilson action (left) and the Mtm-Clover action 
(right). Note that the square mass (a/j,) 2 and the mass (an.) are plotted respectively on 
the x-axis in the two cases. Solid lines are fits to the lattice data. The hatched area 
in the right plot shows the non-perturbative RI-MOM estimate of Zy obtained from 
non-twisted Clover fermions in Ref. [16]. 



a quadratic dependence on the quark mass, due to contributions of 0(a 2 /i 2 ), and 
a linear dependence induced by terms of 0(ci 2 Aqcd A*)- 5 From the results shown 
in Fig. 3, it can be clearly seen that terms of 0(cl 2 Aq C d /-0 are dominant in the 
Mtm-Clover case, while the quadratic contributions of 0(a 2 fi 2 ) are larger in the 
Wilson case, except that at very small values of the quark mass. 

In both the Wilson and Clover cases, we can correct for mass-dependent dis- 
cretization effects by extrapolating the results for Zy to the chiral limit. In this 
way, we obtained the estimates 

z {Wilsan) = a643(1) ) Z (Clover) = QfQg^ (g) 

Our result for 2^^ son ) agrees with the one reported in Ref. [22] , namely 2^ l/son ) = 
0.6424(4). 

The determinations in Eq. (8) can be also compared to the estimates of Zy 
obtained in the corresponding non-twisted theories (for which /i = 0), since renor- 
malization constants are independent of the quark mass. For the standard Clover 
action, non-perturbative determinations of Zy have been obtained using several 

5 We thank R. Frczzotti and G. Rossi for having drawn our attention to this point. 
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approaches. In Ref. [16], by using the RI-MOM method the estimate Zy over = 
0.772(2) (2) has been obtained, in good agreement with our result in Eq. (8) (this 
comparison is also illustrated in Fig. 3). Similar estimates have been obtained 
with the other approaches [23-25]. The determinations of Z v with standard Wil- 
son fermions are affected, instead, by much larger uncertainties, since the action is 
not improved at 0(a). Typically, the non-perturbative estimates lie in the range 
2?(Wiison) ^ 0.6^-0.8 [25-27]. At variance with these estimates, however, we remind 
that our determination in Eq. (8) is improved at 0(a). Another interesting com- 
parison is provided by the estimates Zy Vilson ' > ~ 0.63 and Zy" lover ^ ~ 0.81 obtained 
using one-loop tadpole-improved boosted perturbation theory [19, 20]. 6 



3.2 Light quark masses 

In the quenched approximation, one of the main advantages of tmLQCD is the 
possibility to perform simulations at small values of the quark masses, without 
encountering the problem of exceptional configurations. This can be very ben- 
eficiary for the extrapolations towards the chiral limit. A notable example, in 
this respect, is the calculation of the light quark masses and, in particular, of the 
average up/down quark mass. 

Moreover, the determination of quark masses with Mtm-LQCD is simpler than 
with the standard Wilson and/or Clover actions. To emphasize this point we 
remind the reader that the renormalized quark mass in Mtm-LQCD is obtained 
from the twisted mass parameter fi as rh q = Z^fi, with Z^ = 1/Zp. Thus fi plays 
a role similar to the bare m AWI in the standard, non-twisted, case. But contrary 
to m AWI , which is obtained from the ratio of the correlation functions and thus 
suffers from both statistical and systematic errors, the bare twisted mass fi is 
known without uncertainties. In addition to the automatic (9(a)-improvement of 
the results obtained by using either the Mtm- Wilson or the Mtm-Clover action, 
this makes the calculation of the light quark mass particularly attractive. 

To exemplify the benefits of using the Mtm-LQCD we now apply the same 
strategy explained in detail in Ref. [21], and determine the values of the strange 
and the average up/down quark masses. Concerning the renormalization constant, 
in the Mtm-Clover case we use the non-perturbative estimate Zp I//MOM (2 GeV) = 
0.525 [16], which corresponds to Z]F(2 GeV) = 0.624. In the Mtm- Wilson case, 
instead, we will use the one-loop boosted perturbative result, namely Zp S (2 GeV) = 

6 We choose the boosted coupling obtained by inverting the perturbative series of the logarithm 
of the plaquette, 



ln(P) =-i.9 2 (3.40/a) 



1 - 1.1905 



5 2 (3.40/a) 



47T 



(9) 
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Mtm-LQCD 


m e (MeV) 


m s (MeV) 


Mtm-Clover, m^ r CAC 


4.3 (2) 


105 (5) 


Mtm-Clover, mf™ 


4.2 (2) 


103 (5) 


Mtm- Wilson, m°f 


4.2 (2) 


103 (6) 



Table 2: Values of the average up/down (mi) and of the strange (m s ) quark masses 
renormalized in the MS scheme at the scale of 2 GeV as obtained using the various 
Mtm-LQCD actions at (3 = 6.0. 

0.596. 

The physical values of strange and of the average up/down quark masses have 
been extracted by fitting the results for the pseudoscalar meson masses collected 
in Table 1 as a function of the twisted quark mass. We find that these results are 
consistent with a quadratic dependence of the form 

{aMpf = Pi + P 2 M + P 3 (a/i) 2 • (10) 

The coefficient Pi parameterizes small, though non- vanishing, (9(a 2 )-discretization 
effects: we obtain Pi = {6(2), 5(2), 3(2)} • 10 -3 in simulations a), b) and c) re- 
spectively. In the range of our simulated quark masses we do not observe any 
deviation from the simple quadratic dependence of Mp on the quark masses which 
could signal the contribution of non analytical chiral logs. In Fig. 4 we show 
that the Mtm-LQCD indeed allows a much better control over the extrapolation 
to the physical pion mass than it was the case with the standard non-twisted 
Wilson/Clover action. 

Our final results for the average up/down and strange quark masses are col- 
lected in Table 2. The determinations obtained with the various Mtm-LQCD 
actions are consistent with each other and in agreement with most of the existing 
quenched estimates available in literature (for a recent review see e.g. Ref. [28]). 

3.3 Isospin breaking effects 

The twisted mass formulation of LQCD breaks both parity and flavour symmetries 
at finite lattice spacing. At maximal twist, these breakings are expected to be of 

CV). 

A clear manifestation of isospin breaking effects is a difference between the 
neutral and the charged pion mass. In the Mtm- Wilson case, this splitting has 
been investigated in Ref. [7]. The main observations are: i) isospin breaking effects 
are numerically large, even for reasonably fine lattice spacings; quantitatively, 
it is found 7q (M^ — M^ + ) = c(a/r ) 2 with c ps 10. ii) The contribution of 
disconnected diagrams helps in significantly reducing this splitting; without adding 
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Fig. 4: Pseudoscalar meson masses squared, in physical units, as a function of the 
quark mass m q renormalized in the MS scheme at the scale u = 2 GeV. The dashed 
line represents the result of the fit to Eq. (10). In the Mp values shown in the plot 
the pure lattice artifact contribution represented by the coefficient P\ in Eq. (10) has 
been subtracted. Solid lines indicate the physical values of the pion and kaon masses. 
The results obtained with the non-twisted Clover action in Ref. [21 ] are also shown for 
comparison. 



the disconnected piece, the value of c is about twice bigger, iii) The optimal choice 
of the critical mass, which minimizes parity breaking effects, does not help in 
reducing discretization effects related to isospin breaking. 

We cannot make the study along the lines presented in Ref. [7] , since we explore 
the benefits of using the Mtm-Clover action at single lattice spacing. Furthermore, 
we neglect the contribution of disconnected diagrams, the calculation of which is 
computationally expensive and statistically noisy. One observes that the pure con- 
nected contribution provides by itself the neutral-charged mass splitting between 
pseudoscalar mesons composed by different, though mass-degenerate, isospin dou- 
blets, namely P l = q^^q' with q ^ q'. In the absence of isospin breaking effects 
the neutral and charged mesons would be degenerate, and thus their splitting still 
represents a measure of isospin breaking discretization effect of 0(a 2 ). 

In Fig. 5 we compare the results for the neutral-charged mass splitting obtained 
by using the Mtm- Wilson and the Mtm-Clover actions, with the optimal and the 
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Fig. 5: Relative mass splitting between charged and neutral pseudoscalar mesons as ob- 
tained by using the Mtm-Wilson and the Mtm-Clover actions. The results are presented 
as a function of the bare quark mass. For both actions, the calculation of the neutral 
meson mass only includes the contribution of connected diagrams (see text for explana- 
tions). 



PCAC determination of the critical mass respectively. In the Mtm-Clover case, 
the results obtained with the two different choices for the critical mass (m^ r CAC and 
m^ on ) do not show significant differences. From Fig. 5, it appears that in the Mtm- 
Clover case the isospin breaking effects are significantly reduced with respect to the 
situation in which the Mtm-Wilson action is used, at least for the pure connected 
contribution considered in this study and in the quenched approximation. 



4 Conclusions 

In this paper we have presented an exploratory numerical study of Mtm-LQCD 
with the Clover term in the action. We have performed quenched simulations with 
both Wilson and Clover twisted fermions at /3 — 6.0, with pseudoscalar meson 
masses ranging from about 1.1 GeV down to approximately 280 MeV. We have 
investigated, in particular, the role of the Clover term in reducing the large (IRD) 
cutoff effects which affect, in general, simulations with Mtm-LQCD when the quark 
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mass becomes of 0{aI& CD ). As anticipated in Ref. [8], we find that adding the 
Clover term to the twisted mass quark action is very effective in reducing these 
effects, as much as using Mtm- Wilson fermions with the optimal choice of the 
critical mass. 

By using Mtm-Clover fermions, we have also performed other interesting nu- 
merical studies. We have computed the renormalization constant Zy, finding 
good agreement with the estimates obtained with the standard Clover action. 
We have presented a calculation of the strange and the average up/down quark 
masses, which uses non-perturbative renormalization and much lighter valence 
quark masses that could not have been simulated with standard Wilson/Clover 
fermions, due to the problem of exceptional configurations. Finally, we have cal- 
culated the contribution of connected diagrams to the charged-neutral pion mass 
splitting, finding that the use of the Mtm-Clover action is also beneficial in re- 
ducing cutoff effects related to isospin breaking, at least for this pure connected 
contribution in the quenched approximation. We believe that all these results 
encourage to pursue this project in exploring the unquenched QCD dynamics at 
small quark masses. 
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